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ABSTRACT 

A posteriori anisotropy study of ultra-high energy cosmic rays (UHECRs) with the Pierre Auger Ob- 
servatory (PAO) has shown evidence of excess of cosmic ray particles above 55 EeV within 18° of the 
direction of the radio galaxy Centaurus A. However, the origin of the excess remains elusive. We simulate 
the propagation of different species of particles coming from the direction of Centaurus A in the Galactic 
magnetic fields, and find that only particles of nuclear charge Z < 10 can avoid being deflected outside of 
the 18° window of Centaurus A. On the other hand, considering the increasingly heavy composition of 
UHECRs at the highest energies measured by PAO, a plausible scenario for cosmic rays from the direction 
of Centaurus A can be found if they consist of intermediate-mass nuclei. The chemical composition of 
cosmic rays can be further constrained by lower-energy cosmic rays of the same rigidity. We find that 
cosmic ray acceleration in the lobes of Centaurus A is not favored, while acceleration in the stellar winds 
that are rich in intermediate-mass nuclei, could meet the requirement. This suggests that the observed ex- 
cess may originate from cosmic ray accelerators induced by stellar explosions in the star-forming regions 
of Centaurus A and/or the Centaurus cluster located behind Centaurus A. 

Subject headings: cosmic rays 



Introduction 



Despite nearly one-century of effort, the origin of ultra-high energy cosmic rays (UHECRs) still remains un- 
known. It is believed that the detections of anisotropy in the arrival direction of UHECRs can provide a useful clue for 
recognising the sources. The Pierre Auger Observatory (PAO) has detected 69 comic rays with energy above 55 EeV 
between January 2004 and December 2009. The correlating fraction of the arrival directions which are closer than 
3.1° from the position of an AGN within 75 Mpc is 38% (i.e. 27 out 69 events), while the expected fraction is 21% for 
an isotropic distribution of these cosmic rays. The significance of the anisotropy has decreased compared to previous 



2007 



The Pierre AUGER Collaboration et al 



estimates w ith smaller datasets (The P ierre Auger Collaboration et al 

20081 l20ld) . On the other hand, it is reported that 13 events are located within 18° of Centaurus A (Cen A) out 



of a total of 69 events in the whole sky a bove 55 EeV, whereas 3.2 events would be expected from this region 
dThe Pierre AUGER Collaboration et aDboiol) . implying that the clustering is not a statistical accident. 



It has long been proposed that AGNs could accelerate partic les to ultra-high energies (e.g. lMannheimlll995 



Boldt & Ghosh 



1999[ Berezinskv et al.1 120061 IPermer et al.l 12009b . Particularly, Cen A, as the near est AGN (d 

„t „i Horn n> * ik.,„ — : *i » »;„i f itud/^d„ / — . — n Ji n- 



3.8 Mpc, Harris et a llboiO ), was suggested by various authors as one potential source of UHECRs (e.g. Cavallol 19781 
Romero et all 1 19961: knchordoqui et al.lboOU lHardcastle etaDl2009i: bureev & Troitskylboiol: |Pe'er&Loebll201ll) . 
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However, there are two issues remaining to be solved before concluding that these UHECRs are from nearby AGNs. 
First, most of the ne arby correlating AGNs are subluminous, so t hey may not be powerful eno ugh to accelerate protons 
to 10 20 eV (see e.g. iWaxmanll 1995c |Piranll2010t IWaxmanll201 ll but see iPe'er & LoebluOl lb . requiring intermediate 
mass nuclei or heavy nuclei as the composition of UHECRs in order to reach these energies. The other problem is the 
influence of magnetic field on the trajectories of cosmic rays. Both extragalactic and Galactic magnetic fields (EGMF 
and GMF) can deflect UHECRs, preventing them from pointing back to their birth places when observed from Earth. 
So the apparent arrival direction is not necessarily the real direction of the source and this influence is more serious 
for heavier nuclei. Recent measurements of the maximum air shower elongations < X max > and their fluctuations 
RMS(< X max >) by the PAO indicat e that the cosmic ray spectrum is gradually dominated by intermediate or heavy 



mass nuclei up to E ~ 4 x 10 19 e\l*l A braham et alJboiol) . although this claim depends on the hadronic interaction 



model at high energies which is not well-known at present 0. Interestingly, Hooper & Taylor (2010) find that good 
fits to the < X max > and RMS(< X max >) data can be found for the case in which the sources accelerate primarily 
intermediate mass nuclei. 

In this paper, we examine the effect of Galactic magnetic fields on the deflection of UHECRs of different species 
that originate from the direction of Cen A and study the constraints on the composition of the UHECRs imposed by 
the excess in the direction from Cen A. Since the EGMF is poorly known, its impacts on UHECR deflection is not 
well understood. However, in general it is found that the deflections by th e EGMF tend to be along and within the 



cosmic large scale structure of the galaxy distribution ( Giacinti et al. 201ob . For these reasons, we neglect the effect 



of the EGMF in this work and only consider the effect of the GMF on the deflection of UHECRs. 

The rest of the paper is organized as follows. In §2, we briefly introduce the GMF models employed in this paper, 
as well as our method for simulating the propagation of UHECRs. We show our results and their implications in §3. 
In §4, we further constrain the chemical composition of UHECRs based on our analysis of the lower-energy cosmic 
rays with the same rigidity. We discuss the possible sources of the excess and give our conclusion in §5. Throughout 
the paper, we use eV as the unit of particle energy and use c.g.s units for other quantities and denote by Q x the value 
of the quantity Q in units of 10*, unless specified otherwise. 



2. Simulating propagation of UHECRs in the GMF 



The Galactic magnetic fields are generally described to be composed of a regular component in the dis k and a 
large- scale field in the halo. The previous models for the magnetic field in t he disk are either axissymmetric (e.g. lstanev 
Il997h or bi-symmetric (e.g. lHan et al. | | 1999: Tinvako v & Tkach ev 2002) However, neither of these simple models 
agrees very well with observation s ( Han et all2006l) . Thus, in the present work, we employ t wo up dated models, one 
described in iGiacinti et al.l (1201 C ), which is based on the model built by iProuza & Smfdal d2003l) (hereafter the PS 
model), and another developed in [jiang etalJlhoid) (hereafter the J model) respectively. Both of these models contain 
a disk component and a halo component which consists of a toroidal field and a poloidal field. The halo component 
in these two models are basically the same. The field configurations are based on the anti-symmetric rotation measure 



'Due to lack of statistics, all the events with energy larger than ~ 4 X 10"eV are used in the last bin to get the composition. So we can not get 
any information of X m ax 

and RMS« ^max >) at the highest energy (i.e. < 10 20 eV). In this paper, we assume the tendency which is shown in the 
lower energy range continues to the highest energy. 

2 It should be noted here that observation of High Resolution Fly's Eye Experiment (HiRes) did not show any correlation between arrival direction 
of UHECRs and nearby AGNs, and observation of Telescope Array (TA) can not distinguish a n isotropic distribution from an anisotropic distribution 
of UHECRs arrival directions in current statistics (jAbbasMj^il. 2004 ; Tsunesada^^j20H^. Additionally, both of their elongation measurements 
prefer a pure proton composition in high energy end lATbTs^t^n2^^0Tsune^ad^et^lT2*Olll) . Our paper is based on the observation of the PAO. 
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(RM) s ky revealed by the extragalacti c radio sources dHanetalJI 19971 19991) and the vertical filaments in the Galactic 
center jYusef-Zadeh et all 19841 120 04). while their strengths are selected to meet the features of the obser ved filaments 



in the Galactic center ([Morris & Serabvnl 19961) and the vertical field component in the vicinity of the Sun ( Han & Qiao! 
1994). For the disk component, there is a difference between these two models. In the PS model, the disk component 
is developed from a basic conservative model constructed by Han & Oiaol ( 1994 ). which was based on the Faraday- 
rotation measurement of 134 pulsars, with the as sumption of two logari thmic spiral arms and a constant pitch angle, 
showing bisymmetric (BSS) magnetic field (see lProuza& Srmdall2003l) . While in the J model, the disk component 
is founded upon the RM of pulsars by Han et al. ( 20061). and adopts a four-arm spiral structure, whose configuration 
can be depicted by an Archimedes spira l ( Hou et al.ll2009 ). with the magnetic fields reversed from arms to inter-arms 
( Jiang et allboiO ; Nota & Katgerj2o"lol) . We present the configurations of the disk magnetic field of these two models 
in Figure 1 to show the diff erence. A det ailed description of both the PS model and the J model can be found in 
Giacinti et al. ( 2010 ) and ljiang et al. (2010) respectively. 



Since the deflections of UHECRs caused by the small scale random component at the e nergies we consider 
are smaller or at most com parable to that induced by the regular GMF even for iron nuclei dGiacinti et al. 2010[ 
Tinvakov & Tkachevll2005l) . the random component can be ignored as far as only the largest deflection angle is con- 
cerned in the present study. 



The propagation of cosmic ray particles can be described by 

dp/dt = Zev x B 
dx/dt = v 



(1) 



where e is the charge of the electron and c is the speed of light in the vacuum, p, v and x are the momentum, velocity 
and the spatial coordinates of the particle respectively, and B is the magnetic field. Here we do not consider any 
changes of the parti cle during propagation, such as fragmentation or cooling processes. Following the method used in 
previous works (e.g.lstanevll997[|Harari et al.ll999[ iTinvakov & Tkachevl2002l; IProuza & Smidal2003 ; Giacinti et al. 
20ld Jian g et al we propagate an anti-proton or anti-nucleus from the earth along certain direction, tracking it 



until it reaches the border of the Galaxy, and recording its position and velocity direction at that moment. A proton 
or a nucleus which enters the Galaxy at the same position and with the same direction of its velocity will travel back 
to the Earth along the same path. The spatial extent of cosmic rays emitted by Cen A is assumed to be disk-like, with 
a typical radius of 5° on the celestial sphere. Thus, particles which enter the Galaxy with a direction closer than 5° 
from the center of an object are regarded as coming from that object. 

To study the influence of the GMF on the arrival directions of cosmic rays coming from the direction of Cen A, 
we isotropically propagate and backtrace 10 6 particles with energies from 10 19 75 eV (~ 55 EeV, the threshold energy 
for the clustering) to 10 20 15 eV (~ 142 EeV, the highest energy the PAO ever detected) with an increment of 10 005 
in energy for different species of particles separately. There are two important parameters for each particle. One is 
the angle between the particle's arrival direction and the direction of Cen A, which we denote as 9. The other one 
is the angle between the direction of the particle's velocity when it enters the Galaxy and the direction of Cen A, 
denoted by a. Figure 2 and 3 show the Galactic coordinates of arrival directions of particles coming from Cen A (i.e. 
a < 5°) for the two GMF models respectively. Each dot represents a single particle, and the different colors represent 
their energies. The black filled star is the location of Cen A (I = 309.52°, b = 19.42°) and the black solid curve is the 
projection on the celestial sphere of a circle centered on Cen A with a radius of 18°. The small open circles represent 
the observed events and the size of the circles is proportional to their energies. The 13 magenta circles represent the 
events correlated with Cen A. One can see that, in both GMF models, as expected, the smaller the particle's rigidity 
(i.e. E/Z, where E and Z are the energy and the nuclear charge of the particle respectively), the larger the deviation 
from the original direction. For heavy nuclei, such as iron or calcium nuclei, the particles are severely deflected even 
at the highest energy, while for lighter nuclei such as proton or helium nuclei, the particles can reach the earth without 
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being much deflected. 



3. Effects of GMF on the arrival directions of UHECRs and its implications 

The observed 13 cosmic ray particles from the direction of Cen A spans the energy range from 10 19 75 eV to 
10 19 ,95 eV. The influence of GMF on the arrival directions of cosmic ray particles of different energies can be studied 
quantitatively if their energy spectrum is given. We assume that these particles are accelerated in Cen A to a power-law 
distribution with index s, and keep the same spectrum whilst propagating to the border of the Galaxy, since the distance 
of Cen A to us is much shorter than the attenuation length. Through our simulation of the propagation of these particles 
in the GMF, we get a mock distribution of arrival directions for particles of different species and different energies. 
For convenient comparisons with observations, we define two parameters: one is 0(E), defined as the number ratio of 
particles of species i and energy E that reach the Earth when the deflection by the GMF is considered to those that reach 
the Earth with out considering the GMF deflection. This f actor is also known as the magnetic lensing amplification 
factor (see, e.g. lHarari et al .fcoOOl |2002[ iGiacintiet alboioh . Another parameter is iji(E), defined as the number ratio 



of particles that reach the Earth from a direction within 18° of Cen A (i.e. number of particles with 9 < 18 °) to those 



that r each the Earth from all directions. In the calculations, we consider the relative exposure of the PAO (iSommers 
2001). The corresponding ratios that are integrated over energy from 10 19 75 eV to 10 19,95 eV (approximately the energy 



range of the observed 13 events) can be obtained by 



and 

V ' - JCi(E)§dE ■ (i) 

Table 1 lists the values of these two ratios for different particle species for both the PS model and the J model. From 
the values of Q, one can find that the flux of light particles are barely influenced by the GMF due to their high 
rigidity. Here, the flux of intermediate mass nuclei such as oxygen and silicon are magnified while that of heavy nuclei 
are slightly demagnified. The values of 77, are smaller as the nuclei get heavier, implying that heavy nuclei will be 
deflected severely from their original directions by the GMF. 

Take iron nuclei as an example case. If the excess events from Cen A are dominated by iron nuclei, the all- 
sky event number of iron nuclei contributed by Cen A should be 13/fjp e = 1300 for the PS model, far more than the 
observed event number which is 64 in the same energy range. On the other hand, it has been shown that it is highly 
unlikely for all the observed events with energy > 55 EeV to come solely from a single source, because simulations 
of the arrival directions for such a case indicate that it would di splay an apparently dipolar pattern, which is different 



from the observed distribution even in the most favorable cases ( Giacinti et alfcoi 1 ). Thus, we draw a circle centered 
at Cen A with a radius equal to the largest deflection angle (the pale solid circle in Figure 2 and 3), and regard the 
observed event number in the same energy range that the circle covers as the maximum event number (N mdx j) that 
Cen A can contribute to. For instance, for iron nuclei in the PS model, the maximum event number is 33. In some 
cases (e.g. for light nuclei), the pale circles shrink into the 18° window of Cen A, and we set A^ ma x to be 13 in such 
cases. For iron nuclei in the PS model, they could contribute at most AfcenA.Fe = A^ max .Fe X fjpe = 0.33 events in the 
18° window of Cen A. Thus, as mentioned above, we disfavor iron nuclei as the dominant composition of the excess 
events. 



Similarly, we also calculate the maximum event number that other species of particles can contribute for both 
GMF models. The results are shown in Table 1. One can see that the observed excess events around Cen A can be 
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reconciled to particles with Z < 10 in the both models. Since heavy nuclei such as iron will be deflected severely, 
deviating far from their original direction, they can be ruled out as the main composition of the excess events from 
Cen A. We note that the distribution of simulated particles within the 1 8 window is concentrated in some regions when 
only the regular component of the GMF is considered, which is not as scattered as the observed distribution. Since 
we aim to obtain constraints on the composition of excess cosmic rays, for which only the largest deflection angle is 
concerned, modeling the actual distribution is beyond the present work. However, despite this, we do note that the 
incorporation of random components of the GMF and/or appropriate EGMF would likely was h out the distributio n 
obtained in the present simulation, making it look more like the observed event distribution (e.g. lYiiksel et al.ll2012l) . 



4. Constraints on the source composition with lower-energy cosmic rays 



The shower profile of UHECRs measured by PAO implies that the average chemical composition of cosmic rays 
consists of heavy or intermediate-mass nuclei, most notably at the highest energy bin, as long as the properties of 
hadronic interactions do not change significantly at such high energies. If the measured all-sky-averaged composition 
roughly reflects the average composition of particles from the direction of Cen A, light particles, such as protons and 
helium nuclei, can not be the dominant component of these events. Such a heavy or intermediate-mass nuclei com- 
position could be achieved only if 1) the accelerated material is rich in intermediate-mass or heavy elements; or 2) 
the acceleration ability of the source is limited so that only intermediate-mass or heavy nuclei can be accelerated to 
energies > 10 19 75 EeV (i.e., the source is not powerful enough to accelerate light particles to energies > 10 19 75 EeV). 
In the latter scenario, the composition in the energy range from 10 19 75 eV to 10 19 95 eV is dominated by intermediate- 
mass nuclei or heavy nuclei, regardless of the chemical composition of the accelerated material. However, one should 
be cautious of lower-energy lighter cosmic rays that have the same rigidity, as they will follow the same trajecto- 
ries and could produce stronger anisotropy in the direction of Cen A at lower energies if their abundance is not low 
( Lemoine & Waxmanll2009 : The Pierre Auger Collaboration et alfeoi 1 ). We will show below that the chemical com- 
position of the accelerated material in this scenario must satisfy some constraint as well. 

A number of a cceleration sites of Cen A, such as inner ie ts/lobes ("see. e. g. | Dermer et al . 2009; Rieg er & Aharonian 
2009i lHondall2q09h. the north middle lobe (N ML) (see, e.g. iRomero et al.ll 19961) as well as the giant lobes (see, e.g. 



Hardcastle et alj|2009t lo T Sullivan et alj|2009b . have been sugg ested for UHECRs§ Although there are no precise 
measurements of the element abundance of these possible acceleration regions in Cen A, we can roughly estimate the 
abundance from various observations in these regions. In the inner lobe case, UHECRs could be accelerated by shocks 
driven by the inner lobes inflating into the interstellar medium (ISM). So the chemical composition of the accelerated 
UHECRs in this case should be similar to that of the ISM. Deep Chandra observations of some nearby gas-rich ellip- 
tical galaxies show a nea r-solar metallicity of ISM generally, except an unequivocally sub-solar oxygen abundance 
(< 0.5Z Q , see Kimlboii , and reference therein). For the NML, the spectra of five associated X-ray knots can be well 
fitted by sub-solar element abundance ( Kraft et al]l2009 ). implying a sub-solar metallicity for NML. As for the giant 
lobes, since they extend out to 600 kpc from the center, the ambient medium should be intergalactic medium (IGM), so 
the chemical composition of accelerated UHECRs is probably identical to that of IGM, whi ch is found to be sub-solar 
metallicity both observationally and theoretically (e.g. lRichter et al J2009t Barai etal Eoi lj). Here we assume the solar 
ab undance composit i on for cosmic rays accelerated in these regions and adopt the solar element abundance reported 
in lLodders & Palmd (120091) in the following calculation. For simplicity, we attribute elements with 6 < Z < 10 to 



3 The central black hole with a strong magnetic field has also been suggested as a possible accelerator (e.g. iNeronov et al' though other 

works have highlighted th at the maximum energy may be somewhat limited, even for the of case of nuclei, see e.g. jLemoine & Waxmarjlaxwl 
iRieger & Aharonia r]|2009l) . For this reason we do not consider it in the present work. 
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oxygen, with 11 < Z < 16 to silicon, with 17 < Z < 23 to calcium, and with 23 < Z < 26 to iron, and neglect heavier 
elements. Then the relative mass abundance is M H : M He : M '■ M si : M Ca : M Fe « 60 : 24 : 1 : 0.16 : 0.014 : 0.12. 

Since cosmic rays of solar abundance composition are dominated by light elements such as hydrogen and helium, 
the condition 1) can not be satisfied and hence the condition 2) should be taken into consideration. We assume that the 
maximum and minimum energies of particles for different species are rigidity-dependent (i.e. £ ma x,£min °c Z), and the 
spectrum for nuclei with nuclear charge Z, is described by 

dN i /dE = f i N (E/E )- s E <E<E 

f.max ■ (4) 

where /* is the relative number abundance of element i at a given energy E, and iVo and Eq are used for normalization. 
Above E m - dx or below E m i n , an abrupt cutoff in the spectrum is assumed. Thus, the relative number abundance /; relates 
to the mass abundance by 

Mi : Mj = ^ A^dE : ^ Aj^dE =A i Z}~ s f i : AjZ^fj (ij = H,He,C, ...)■ (5) 

For a mixed composition, the maximum number of events within 18° of Cen A is given by 

Efw^Q(EME)§dE 

^VcenA.mix = /.10"- 95 eV ~JZ * Mnax.mix (6) 

Sil0 19 - 75 eV 0(E)-j£dE 
i 

where N mdx . m i x is the maximum event number that could be contributed by Cen A, as defined in §4, for such a mixed 
composition. The value of N maXtmix can be taken as 30, since one can find from Figure 2 and Figure 3 that for most 
species this value is roughly < 30 (see discussion in Section 3). Similarly, for a mixed composition, the average atomic 
mass of cosmic ray events contributed by Cen A is 

v r 10flM ° y A-e-(E)&dE 

i 

One can see that both A^enA,mix and < A > are related with the maximum acceleration energy E i max through the 
spectrum dN,(E)/dE. As the maximum energy is rigidity-dependent, E j mm = ZjE p max , where E p max is the maximum 
accelerating energy for protons. Any appropriate value of E pAWdx should lead to A^enA.mix ^ 10 (assuming 3.2 out 
of 13 events are from isotropically distributed sources in this region). On the other hand, the fits to both < X max > 
and RMS< X mdx > have shown that the all-sky-averaged composition of amvingUHECRs may be nitrogen-like 
(ITaylorll201 II) . if only a small spread in composition exists at energies ~ 10 19 5 e\o Under the assumption that the 



events contributed by Cen A have the same composition as the measured all-sky composition, we conservatively set 
< A >> 10 to be another necessary condition. 

Figure 4 shows the maximum number A^cenA.mix an d the average atomic mass < A > of cosmic rays from the 
direction of Cen A as a function of E p max for solar abundance composition of cosmic rays. The red lines represent the 
results for the PS model of the GMF, while the blue ones represent the results for the J model. One can see that in both 
GMF models, low E pfam (e.g. < 10 19 eV) can be excluded, because in this case only heavy nuclei can be accelerated up 



More generally speaking, it is difficult to explain < X mllx > and RMS(< X max >) simultaneously by any composition because of the low value 
of the measured RMS(< X max >). Indeed, a larger RMS(< X max >) is predi cted by transition models in which the composition changes fro m a light 
to a heavier composition, as is suggested by the trend in < X max > (see e.g. lThe Pierre Auger Collaboration et al.fetH lUuirich e t al. 2009). 
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to the required energy, while for heavy nuclei, due to their low-rigidity, most of them are deflected into other directions 
or even can not reach the Earth. So E pm!lx should be large enough to allow at least intermediate-mass nuclei to be able 
to reach 10 19,75 eV. On the other hand, because of the large abundance of hydrogen or helium elements relative to 
oxygen (and heavier elements), the average chemical composition will become too light once FHe.max > 10 19,75 eV. As 
a consequence, we have the pale shaded region and the dark shaded region (including the pale one) representing the 
range of E pMi - dx in which both requirements can be satisfied for the PS model and the J model respectively. Note that 
the limit of E pfasK obtained here is consistent with the theoretical acceleration limit of Cen A (Lemoine & Waxman 
2009; Piran 2010), although the acceleration ability of Cen A is still under much debate (Rieger & Aharonian 2009; 
Pe'er&Loeb 2011). 

Lower-energy cosmic rays that have the same rigidity (E/Z) will travel along the same paths. The chemical 
composition of cosmic rays can be constrained by requiring that they do not produce excess in cosmic-ray flux in the 
direction of Cen A and in the all-sky flux at corresponding energies. As shown above, in the viable range of E p . m:ix 
obtained, the observed excess events from the direction of Cen A are very likely to be dominated by intermediate 
mass nuclei with 6 < Z < 10, so we take oxygen nuclei as the dominant particles in the following calculation. Given 
that the ratio between the number of excess cosmic rays from the direction of Cen A and that of the whole sky is 
10: 64 in the energy range from 10 19 75 eV to 10 19 95 eV, cosmic rays originating from Cen A should account for a 
flux F ~ 10/64F ob (10 19 - 75 - 19 95 eV) within the 18° window around Cen A and an all-sky flux of F . a ii = F /r] in 



sxpfiQiMs ev ^JeT^E i s me measured all-sky flux, dN° b jdE is 



the same energy range, where /r°b(io 19 ' 75 " 19 ' 95 e y) = A 
the measured differential flux, and A exp is the total exposure of the PAO. So particles of species i that have the same 
rigidity in the energy range from JMO 19 75 eV to JMO 19 ' 95 eV should produce a flux of F = ^ (J^) 1 ~ s Fo inside the 18° 
window of Cen A and an all-sky flux of F !;a u = ^-( J^) 1_s Fo. a ii, respectively. 

Since the search at lower energies has not revealed any significant anisotropy signals around the Cen A region 
dThe Pierre Auger Collaboration et al.ll2.Ql ID . the flux F; of lower-energy cosmic rays should not be higher than the 
expected background flux, i.e. 



i -j 



Id 1 



A) 



19.75 



10 19 - 95 ) eV 



(8) 



where x ~ 0.0466 is the fraction of the exposure of the PAO within the 18° window of Cen A (IThe Pierre Auger Collaboration et al 



201 lb . It should be pointed out here that the above inequality is a conservative constraint, since the anisotropy signal 
can be fairly strong when the source counts are comparable to the background counts, especially if the number of 
background counts is large. Note that this method is only applicable to lower-energy, lighter cosmic rays. 

One can also get constraints on the cosmic ray composition by requiring that lower-energy all-sky flux produced 
by sources in the direction of Cen A should be lower than the measured flux in the same energy range, i.e. 



Si (Zi 



So \Zo 



Fo/vo<F ob 



Zo 



(10 



19.75 



lO 19 ' 95 ) eV 



(9) 



The m easured cosmic ray flux in the relevant range can be fitted by a broken power-law, given by ([Abraham et al 
201 Oh 



dN, 



ob 



,E<E h 



(10) 



< IE I M. I - 1 ,E>E h 
where Et, = 10 19 46 eV is the break energy, N\, is the normalized factor, p\ = 2.6 and pi = 4.3 are the power-law indexes 



-8- 



of the two segment, respectively. Then from inequality © we obtain 



fi < fo x 6.4* ( §- 



-* ( l-P2 S 
Z Q ) \\-p X/ 

and from inequality © we obtain 

6.4*, (g)** 

6A *° (I) 



for lighter elements, 



(ID 



for heavier elements, 



££' '' 2 for lighter elements. 



(12) 



Here E\ = 10 eV and E2 = 10 eV. One can find that for lighter elements, which method gives a stronger 
constraint depends only on the value of x and 770. In the case discussed in the present work, x <C r?o, so the constraints 
from the 18° window anisotropy are stronger. But for heavier elements, useful constraints only come from the latter 
method. Assuming s = 2, inequality (11) and (12) give f H < 4.4/ , / He < 2.9/ , f S i < 0.72/ o , / Ca < 0.3 1/ and 
fee S O-17/o, corresponding to M H < 2.2M , M He < 2.9M , M si < 0.72M o , M Ca < 0.3 1M , M Fe < 0.18M o for the 
PS model; and f H < 4.4/ , / He < 2.9/ , / si < 1 .2/ , / Ca < 0.5 l/o and / Fe < 0.29/ o , corresponding to M H < 2.2M , 
Mm < 2.9M , M S i < 1.2M , M Ca < 0.5 1M , M Fe < 0.3 1M for the J model. A larger power-law index s will put 
a more stringent constraint on the relative abundance of lighter particles. One can find that, in both GMF models, 
constraints on heavier elements are easy to satisfy, i.e. the excess in flux at higher energy induced by heavier nuclei with 
the same rigidity can be avoided. But to avert a low energy excess in the Cen A direction made by protons and helium 
nuclei with the same rigidity, a super-solar metallicity especially a super-solar intermediate-mass element abundance 
is re quired. This resu lt is consistent with the analysis from fitting X max and RMS(< X max >) that measured by the 
PAO dShaham & Piranboij . Such a chemical composition is not favored by the candidate cosmic-ray accelerators in 
the radio galaxy Cen A. In other words, if we attribute the 10 excess events to Cen A, a stronger excess would also 
occur at 1-10 EeV. We note, however, that if the composition of these 10 clustering events turns out to be proton or 
helium nuclei, different from the meas ured all-sky- averaged composition, then the above conclusion does not apply 
(for light nuclei composition, see e.g. FargionlEoi 1 , and reference therein). 



5. Discussion and Conclusion 



The required chemical composition of the excess cosmic -rays in the direction of Cen A seems to disfavor the 
above mentioned candidate sites in Cen A. The requirement, however, can be more easily satisfied by UHECR acceler- 
ators induced by stellar explosions since either the stellar wind of massive stars or the exploded ejecta can be rich in in- 
termediate mass nuclei. The chemical co mposition of WR stellar winds is Mu e ■ Mq : Mo : Mx = 0.32 : 0.39 : 0.25 : 0.04 
( Bieg i"ngfl990l : Ivan der Hucht et alii 19861) where Mx is the sum of elements heavier than oxygen, which we treat sim- 
ply as silicon in the following calculation. The composition of the hyper nova ejecta is also rich in oxygen nuclei. The 
numerical modeling of the early spectra and light curve of SN 1998bw dNakamura et al. 2001 ) yields a composition 
of M c : M : M Ne : M Mg 



M Si : M s : M Ca : M Fe = 0.006 : 0.71 : 0.037 : 0.034 : 0.083 : 0.041 : 0.007 : 0.09. We show 



the dependence of Nc en A and < A > on E p msx in the WR stellar wind and hypernova ejecta composition scenarios 
respectively in Figure 5. One can see that as long as E p nax is larger than 10 19 eV, which can be reached in the GRB 
and hypernova scenarios dWang et al. 2007 20081) . both AfcenA and < A > meet the requirements. It has been also 
shown that the spectrum and composition o f UHECRs accelera ted in the WR stellar wind or in the hypernova ejecta 
are compatible with the PAO's observations ( Liu & WangH2012 ). 



Both (long) GRBs and hypernovae are generally believed to trace star formation. Recent observations by the 
Hubble Space Telescope revealed triggered star formation occurred < 10 Myr ago in the inner filament of Cen A 
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dCrockett et al.ll2.012D . although Cen A belongs to elliptical galaxies, which usually have low star formation rates. 
Possible alternative sites are GRBs and hypernovae occurring in the Centaurus cluster, which lies behind Cen A at 
a distance of ~ 50 Mpc. Since there are hundreds of galaxies in the Centaurus cluster, the star formation rate is 
much higher than that in Cen A, and as a result, the rate of GRB or hypernova is much higher. For a distance of 
~ 50 Mpc, intermediate-mass nuclei such as oxygen of > 60 EeV will suffer from photo-disintegration during the 
propagation, so the chemical compositi on at the border o f the Galaxy after propagation may be different from that at 



the sources (also see the discussion in lTaylor et al.ll201 ll) . Figure 6 shows the propagated spectrum of oxygen nuclei 
from the Centaurus cluster. The integrated flux from 10 19 75 eV to 10 19,95 eV are about a factor of 0.15 of the initial 
one. Heavier nuclei at higher energies also suffer from severe attenuation, while helium nuclei and protons at lower- 
energies (~ 10 18 eV- 10 19 eV) are almost not affected. Taking into account such attenuation, therefore, leads to even 
stronger constraints on the element abundance of light particles at the sources, i.e. approximately reducing to 0.15 
times the value obtained in the previous section, leading to the new constraint of Mh < 0.33Mo and Mu e < 0.43Mo- 
Even though this new constraint is harder to satisfy, one finds that the chemical compositions of WR stellar winds and 
hypernovae ejecta can still satisfy such a strict requirement. 

In this paper, we study the origin of the observed excess of UHECRs in the direction of Cen A. First, by simulating 
the propagation of cosmic rays in the GMF, we find that the excess events can not be mainly composed of heavy nuclei 
like iron nuclei coming from the direction of Cen A , because the GMF has a significant influence on their trajectories 
and cause their apparent arrival direction to severely deviate from their original directions. Also, the excess events are 
not likely to be dominated by light particles such as helium nuclei or protons as the measurements of the elongation 
rate of air showers by PAO suggest a heavy or intermediate mass composition. We show that intermediate mass nuclei 
with nuclear charge 6 < Z < 10 are good candidates of the main composition of the excess events. We also show 
that the composition of cosmic rays can be further constrained when the anisotropy of low-energy cosmic rays is 
considered. In order not to produce a significant anisotropy in the direction of Cen A at lower energies, which is not 
observed by the PAO, a low proton and helium abundance (e.g. a supersolar metallicity) is required. None of the 
proposed candidate acceleration sites in Cen A is favored by this constraint. We find that the cosmic-ray accelerators 
arising from stellar explosions such as GRBs or hypernovae are more favorable sources of the excess events because 
high abundance of intermediate mass elements could be possible in these accelerators. Since the event rates of GRBs 
or hypernovae trace the star formation rate, the star formation region in the inner filament of Cen A or in the Centaurus 
Cluster may be the origin of the excess events. However, as a word of caution, it should be remembered that the excess 
is significant only a posteriori. Furthermore we note that due to uncertainties in our knowledge of the hadronic model 
describing cosmic ray atmospheric showers, the GMFs, the chemical environment of Cen A and even the acceleration 
mechanism of UHECRs, the conclusions made in this paper should be regarded as model-dependent. 
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Table 1 : Some useful parameters obtained in our simulation. ( is the averaged number ratio of particles that originate 
from the direction of Cen A and reach the Earth when the GMF deflection is considered to those that reach the 
Earth without considering the GMF deflection, fj is the averaged number ratio of particles that reach the Earth from 
a direction within 18° of Cen A to those that reach the Earth from all directions, and A^cenA is the maximum event 
number of cosmic rays that are confined within 18° of Cen A out of 64 UHECR events. 



Elements 


c 


PS model 

fj 


AfcenA 


C 


J model 

fj 


AfcenA 


H 


1.00 


0.96 


27.8 


0.95 


1 


13 


He 


1.03 


0.93 


27.0 


0.95 


1 


13 





2.12 


0.41 


11.9 


1.84 


0.66 


20.5 


Si 


1.80 


0.02 


0.58 


2.55 


0.25 


8.00 


Ca 


0.86 


0.008 


0.23 


0.79 


0.07 


2.31 


Fe 


0.88 


0.01 


0.33 


0.44 


0.004 


0.13 
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Fig. 1 . — Configuration of the disk magnetic field in the PS model (upper panel) and the J model (bottom panel), 
respectively. The direction of arrows represents the direction of the field while the length of arrows represents the 
strength of the field. In particular, magnetic field > 3/iG is colored as red, > 2/iG yellow, > lfiG green, < 1/iG blue. 
The open circle shows the location of the sun. For clarity, the field at the Galactic center is not shown. The shaded 
region in the bottom panel outlines four spiral arms. 
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Fig. 2. — Arrival directions (in Galactic coordinates) of cosmic rays originating from the direction of Cen A in the 
PS model. A total of 10 6 particles are used in these simulations. One dot represents one particle, and different 
color dots/bands represent different energies (from 10 19 75 eV to 10 20 15 eV with an increasement of 10 ). The black 
filled star is the location of Cen A (/ = 309.52°, b = 19.42°) and the black solid curve is the projection on the celestial 
sphere of a circle centered at Cen A with a radius of 18°. The pale solid curve is the projection of a circle centered at 
Cen A with a radius equal to the largest deflection angle of arriving cosmic rays. The small open circles represent the 
observed events and the size of the circle is proportional to the energy. The 13 magenta ones are those events that are 
correlated with Cen A. See text for more discussions. 
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Fig. 3. — The same as Figure 2 but for the case of the J model. See text for more discussions. 
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Fig. 4. — The maximum number of cosmic rays N cen A from the direction of Cen A and the average atom mass < A > 
versus different maximum proton acceleration energy E p ^ max for solar abundance composition of cosmic rays. The 
average atom mass < A > is shown by the dashed lines with the value indicated by the left-handed vertical axis, while 
AcenA is shown by the solid lines with the value indicated by the right-handed vertical axis. The corresponding values 
of £p, ma x where "Si", "O" and etc are indicated are the lower limits of E p max for which nuclei of such elements can 
be accelerated to 10 1975 eV. The shaded region represents the appropriate range of E Ptmax that satisfy both N cen A ^ 10 
and < A >> 10. The verti cal dash-dotted line shows an approximate theoret i cal acceleration upper limit for proton s 
in Cen A (~ 10 19 eV, e.g. lHardcastle et"ai]|2009t Eemoine & Waxmanll2009l: |pirar]|20ld: bureev & Troitskvlboiol) . 
however, which is still under much debate. See text for more discussion. 
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Fig. 5. — The same as Figure 4 but for different cosmic-ray composition. The top panel corresponds to a WR stel- 
lar wind composition and the bottom panel corresponds to the composition of hypernova ejecta. See text for more 
discussions. 



-19- 




Fig. 6. — The propagated spectrum of oxygen nuclei after traveling a distance of 50 Mpc. The initial spectrum is set 
to be a power-law spectrum with index of -2 and with an abrupt cutoff atl0 21 eV. 



